Draft version January 20, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



DISCOVERY OF A MULTIPLY LENSED SUBMILLIMETER GALAXY IN EARLY HERMES 

HERSCHEL/SPIKE DATA 

A. Conley\ A. CooRAY^'^, J.D. ViEiRA^, E.A. Gonzalez Solares*, S. Kim^, J.E. Aguirre\ A. Amblard^, R. Auld", 
A.J. Baker', A. Beelen^, A. Blain^, R. Blundell^, J. Bock^'^", CM. Bradford^"', C. Bridge^, D. Brisbin^\ 

D. Burgarella^^, J.M. Carpenter'^ P. Chanial^^, E. Chapin^^, N. Christopher^'', D.L. Clements^*", P. Cox^"", 

S.G. Djorgovski^, CD. Dowell^"', S. Eales", L. Earle^ T.P. Ellsworth-Bowers^**, D. Farrah^", 
A. Franceschini^", D. Frayer", H. Fu^', R. Gavazzi^% J. GLENN^'^^ M. GRIFFIN^ M.A. GuRWELL^, M. Halpern", 

E. Ibar,^^, R.J. IvisoN^'''^^ M. Jarvis^"', J. Kamenetzky^*, M. Krips", L. Levenson'''^" R. Lupu"', A. Mahabal^, 

P.D. Maloney\ C. Maraston^", L. Marchetti^'^, G. Marsden^^, H. Matsuhara^^ A.M.J. Mortier^'"', 

E. MURPHY^'^**, B.J. NAYLOr"\ R. NeRI^^ H.T. NGUYEN"'■^ S.J. OLIVER^^ A. OmONT^S M.J. PAGE2^ 

A. Papageorgiou" , CP. Pearson''"''\ I. Perez-Fournon^^"*^, M. Pohlen^, N. Rangwala\ J.I. Rawlings^®, 
G. Raymond**, D. Riechers^-^^, G. Rodighiero'", I.G. Roseboom^'', M. Rowan-Robinson"', B. Schulz^'^*, 
Douglas Scott", K. Scott", P. Serra^, N. Seymour^", D.L. Shupe^'^*, A.J. Smith"-", M. Symeonidis^^, 

K.E. TUGWELL^^ M. VACCARI^", E. VALIANTe", I. VALTCHANOV^'^ a. VeRMA^^ M.P. VIERO^ L. VlGROUX^^ L. WANG^^ 

D. Wiebe", G. WRIGHT^^ CK. Xu^'^*, G. Zeimann^'', M. Zemcov^'^", and J. Zmuidzinas^'^" 

Draft version January 20, 2013 

ABSTRACT 

We report the discovery of a bright (/ (250/ini) > 400 mJy), muhiply-lensed submihimeter galaxy 
HERMES Ji q5751.1-H573027 in Herschel/SPIRE Science Demonstration Phase data from the Her- 
MES project (jOUver et al.ll2010l) . Interferometric 880/im Submihimeter Array observations resolve at 
least four images with a large separation of ^9". A high- resolution adaptive optics Kp image with 
Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of 2; = 2.9575, 
and the lensing model gives a total magnification of /i ~ 11 ± 1. The large image separation allows us 
to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by 
the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody 
fit with an unusually warm dust temperature, 88 ± 3 K. We derive a lensing-corrected total IR lumi- 
nosity of (1.43 ± 0.09) X 10^'^ L0, implying a star formation rate of ~ 2500 Mq yr~^. However, models 
primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full 
optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered 
in early Herschel data, and many more are expected as additional data are collected. 
Subject headings: galaxies: high-redshift — galaxies: starburst — gravitational lensing: strong — 
submihimeter: galaxies 



I alex ander.conley@colorado.edu | 

'"Center for Astrophysics and Space Astronomy, 389 UCB, 
University of Colorado, Boulder, CO 80309 

^ Dept. of Physics & Astronomy, University of California. 
Irvine, CA 92697 

^ California Institute of Technology, 1200 E. California Blvd., 
Pasadena, CA 91125 

* Institute of Astronomy, University of Cambridge, Madingley 
Road, Cambridge CB3 OHA, UK 

^ Department of Physics and Astronomy, University of Penn- 
sylvania, Philadelphia, PA 19104 

® Cardiff School of Physics and Astronomy, Cardiff University, 
Queens Buildings, The Parade, Cardiff CF24 3AA, UK 

Department of Physics and Astronomy, Rutgers, The State 
University of New Jersey, 136 Frelinghuysen Rd, Piscataway, NJ 
08854 

* Institut d'Astrophysiquc Spatiale (IAS), batiment 121, Uni- 
versite Paris-Sud 11 and CNRS (UMR 8617), 91405 Orsay, 
France 

® Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 

^'^ Jet Propulsion Laboratory, 4800 Oak Grove Drive, 
Pasadena, CA 91109 

Space Science Building, Cornell University, Ithaca, NY, 
14853-6801 

^'^ Laboratoire d'Astrophysiquc de Marseille, OAMP, Univer- 
site Aix-marseille, CNRS, 38 rue Frederic Joliot-Curie, 13388 
Marseille cedex 13, France 

13 Laboratoire AIM-Paris-Saclay, CEA/DSM/Irfu - CNRS - 
Universite Paris Diderot, CE-Saclay, pt courrier 131, F-91191 
Gif-sur-Yvette, France 



^'^ Department of Physics & Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver, BC V6T IZl, 
Canada 

Department of Astrophysics, Denys Wilkinson Building, 
University of Oxford, Keble Road, Oxford OXl 3RH, UK 

Astrophysics Group, Imperial College London, Blackett 
Laboratory, Prince Consort Road, London SW7 2AZ, UK 

^'^ Institut de RadioAstronomie Millimetrique, 300 Rue de la 
Piscine, Domaine Univcrsitairc, 38406 Saint Martin d'Heres, 
France 

Dept. of Astrophysical and Planetary Sciences, CASA 389- 
UCB, University of Colorado, Boulder, CO 80309 

Astronomy Centre, Dept. of Physics & Astronomy, Univer- 
sity of Sussex, Brighton BNl 9QH, UK 

Dipartimento di Astronomia, Universita di Padova, vicolo 
Osservatorio, 3, 35122 Padova, Italy 

21 NRAO, PO Box 2, Green Bank, WV 24944 

22 Institut d'Astrophysiquc de Paris, UMR 7095, CNRS, 
UPMC Univ. Paris 06, 98bis boulevard Arago, F-75014 Paris, 
France 

23 UK Astronomy Technology Centre, Royal Observatory, 
Blackford Hill, Edinburgh EH9 3HJ, UK 

2^ Institute for Astronomy, University of Edinburgh, Royal 
Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

25 Centre for Astrophysics Research, University of Hertford- 
shire, College Lane, Hatfield, Hertfordshire ALIO 9AB, UK 

26 Institute of Cosmology and Gravitation, University 
of Portsmouth, Dennis Sciama Building, Burnaby Road, 
Portsmouth POl 3FX, UK 

2^" Institute for Space and Astronautical Science, Japan 



2 



Conley et al. 



1. INTRODUCTION 

The discovery of a population of high-redshift, 
prodigiously star-forming galaxies at sub-millimeter 
wavelengths has revolutionized our understand i ng of 
cosmological star fo rmation (e.g., iBlain et all 120021 : 
iChapman et al.l I2005D . These submillimeter galaxies 
(SMGs) are frequently faint at optical wavelengths due 
to significant extinction, but some have far-infrared lu- 
minosities in excess of 10^'^ L©, and are forming stars at 
> 1000 M© yr~^. They are believed to be the progeni- 
tors o f nearby massive elliptical galaxies (jSwinbank et al.l 
|2008() . yet many of their properties remain mysterious. 

Dusty, star forming galaxies are responsible for 
most of the cosmic infrar ed background (GIB, e.g., 
IMarsden et al. II2009I: iGleim et al.ll201(]D . which contains 
as much en ergy as all of the optical light ever emitted 
by galaxies ()Puget et al.lll996l ). Modelers have had some 
success in fitting the spectral energy distributions (SEDs) 
of SMGs and using this to infer their properties (e.g., 
IRowan- Robinson et al.l [20081 ) . but it is difBcult to study 
the sources that produce the CIB in detail because they 
are individually faint. These efforts are biased towards 
extremely luminous, red galaxies by selection effects, so 
it is interesting to test how well such models describe the 
less luminous Herschel sources selected at shorter wave- 
lengths. 

Confusion noise generally sets the flux limit at which 
individual Herschel sources can be studied. Gravita- 
tional lensing allows this limit to be circumvented. Due 
to t he rapidly rising s ource counts at faint flux densi- 
ties (jGlenn et al.ll2010l) and the negative X-correction at 
sub-millimeter wavelengths, strong lensing is expected 
to be relatively common for SMGs (JBlain 1996). In- 
deed, follow up of early SPIRE data has shown that 
a large fractio n of the bright est sources are lens ed by 
other galaxies (jNegrello et al.ll2010D . iVieira et all H2010f ) 
discovered a population of bright galaxies at 1.4 and 
2mm, and suggested that these are lensed. Lensing al- 
lows us to study the properties of intrinsically fainter 
SMG s at a level of detail that is currently diffic ult other- 
wise ([Swi^inkiriniloil [ivioniiri^ Galaxy- 
galaxy lensing is expected to dominate, with generally 
small image separations, so emission and absorption as- 
sociated with the foreground lens may obscure the SMG 
in the opt ical and near-IR; this is the case for all of the 
sources m iNegrello et aD ([2010| ). 

Here we report the discovery of an SMG system (HER- 
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MES J105751. 1-1-573027, hereafter HLSW-01) multiply 
lensed by a group of galaxies at z = 2.9575 ± 0.0001 
in Science Demonstration Phase Herschel/SPIRE ob- 
servations of the Lockman-SWIRE field as part of the 
Herschel Multi-tiered Extragalactic Survey (HerMES; 
S. Oliver et al. 2011, in preparation), with coordinates 
a = 10h57m51s 6 = 57°30'27" (J2000). A number of 
additional lensed systems are already known in HerMES 
data. The large separation between the images allows us 
to measure the SED of this object across a long wave- 
length baseline. In this letter, we model the optical-to- 
millimeter SED of this object. The lensing model for this 
system, based on high-res olution optical and ne ar-IR ob- 
servations, is presented in lGavazzi et al.l (|2011l hereafter 
Gil). We have also obta ined high-resolution CO line 
maps (jRiechers et al.l[2010l hereafter Rll), and used the 
GO lin e strength distribution to model the molecular gas 
(jScott et al..,2011l . hereafter Sll). 

2. OBSERVATIONS 

HLSW-01 was discovered using observations with the 
Spectral and Photometric Imaging Receiv er (SPIRE, 
i Griffin et al.l [2010) on-board Herschel (Pil bratt et al.l 
I2OIOD . It was selected for further follow-up with Z-Spec, 
a millimeter-band grating s pectrograph at t he Caltech 
Submillimeter Observatory (jEarle et al.ll2006h . based on 
its brightness and blue color (/ (SOOyitm) < /(300/xm)); 
the latter was intended to avoid z > A sources where 
a redshift would be difficult to obtain. Z-Spec gives a 
secure redshift of z = 2.958 ± 0.007 (Sll). Rll use addi- 
tional GO lines measured with the Plateau de Bure Inter- 
ferometer (PdBI), the Combined Array for Research in 
Millimeter- wave Astronomy (CARMA), and Zpcctromc- 
ter on the Green Bank Telescope to refine the redshift, 
yielding z = 2.9575 ± 0.0001. 

HLSW-01 is unresolved in the diffraction-limited 
SPIRE observations (EWnMsso^m = 18.6"). The ex- 
treme brightness (/ (250/zm) ~ 400 mJy) and the mor- 
phology in the optical and near-IR suggested a lensed 
source. A Subaru i image from the SERVS surve£!| 
shows clear evidence of lensing. The source is visible in 
previously obtained gr WHT observations and in data 
from the Spitzer/ SWIRE survey (Surace et al. 2011, in 
preparation) . We obtained observations with the Submil- 
limeter Array (SMA) in compact configuration at 880/im 
(beam FWHM 2.3"), resolving the source into at least 
four components matching the position of several optical 
sources and surrounding a foreground elliptical galaxy. 
We further obtained a Kp adaptive optics (AO) observa- 
tion using NIRG2 on the Keckll telescope and the laser 
guide-star system (jWizinowich et al.l 120061 ) . PdBI CO 
maps, presented in Rll, detect at least the two brightest 
sources and show that they are at the same redshift. The 
photometric re dshift of the central elliptical is 0.60 ±0.04 
(jOvaizu et al.|[200 8i). An optical spectrum obtained with 
the double spectrograph on the Hale telescope has ab- 
sorption features consistent with the photo-z. 

The numbering scheme used to identify sources in this 
letter is shown in Figure[TJ along with as a montage of ob- 
servations at multiple wavelengths. The high-resolution 
Kp image was used to construct a lensing model, which 
was compared with the i and PdBI observations to con- 

3'^ |http: //www, cv.nrao. edu/-mlacy/servs .html| 
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strain differential magnifieation; see Gil for details. This 
model has five lensing masses at the locations of fore- 
ground galaxies, and gives a total magnification factor of 
= 10.86 ± 0.68 for all five detected images based on a 
cored-isothermal model. The velocity dispersion strongly 
suggests that the lenses reside in a massive group of 
galaxies. 

The photometry is summarized in Table [T] At optical 
wavelengths the individual images are blended and par- 
tially resolved, so neither point-spread-function (PSF) 
nor aperture photometry is entirely satisfactory. We use 
aperture photometry with a relatively small aperture ra- 
dius to minimize blending, and compute the aperture 
corrections for the partially resolved sources by convolv- 
ing the lens model to the matching resolution in each 
band; the adjustment to the aperture correction com- 
pared with isolated point sources is only a few percent. 
We exclude photometry of image 1 because it is contam- 
inated by a foreground object, and again use the lensing 
model to correct for the omitted light, which is a ~ 15% 
correction. 

In the 5piteer/SWIRE data, the individual images are 
separable in the Infrared Array Camera (IRAC) bands 
(3.6 to 8/im), and the SWIRE PSF photometry is ad- 
equate for our purposes. We use only photometry of 
sources 2 and 3, again using the lensing model to correct 
for the omitted sources 1 and 4. HLSW-Ol is not present 
in the SWIRE catalog at 3.6/im. The FWHM of the 70 
and 160/im MIPS observations are much larger than the 
separation between individual images, so the catalog flux 
measurement already includes all the images. At 24/im 
the SWIRE aperture is just smaller than the separation, 
so we re-measured the photometry using a larger aper- 
ture. 

The SPIRE fluxes used here are from the HerMES 
SCATvS.l catalog (A. J. Smith et al. 2011, in prepa- 
ration). For the SMA data we extract photometry and 
positions using the CASAF^ imf it task. For Z-Spec, we 
bin the spectrum into 5 100/im bins after masking noisy 
channels and the detected CO l ines. The 20c m photom- 
etry is from the FIRST survey (|Becker et aLl flQQS). 

We place upper limits on the potential foreground con- 
tami nation by scaling t h e observed el l iptica l galaxy SEDs 
froml Pale et all ()2007| ): iTemL et al.l ()2007l ) to match the 
optical magnitudes of the central lensing elliptical, taking 
the highest resulting fluxes in each band as our contami- 
nation limit. This is not relevant at shorter wavelengths 
where the sources are clearly resolved. The potential 
contamination peaks at 6 mJy at 160/im, and 3 mJy at 
70/im. Given the SED, only the potential 70/im contam- 
ination is significant, but we adopt the contamination 
limits as an additional correlated uncertainty at all wave- 
lengths. We take calibration errors, which also affect the 
photometry in a correlated fashion, from the instrument 
documentation. 

3. SED FITTING AND SOURCE PROPERTIES 

We first analyze the SED by fitting simple modified 
blackbody models to the long- wavelength data. It is pos- 
sible that the far-IR emission has a different spatial dis- 
tribution than the shorter wavelength data used to derive 
the lensing model, but the current SMA observations are 

http://casa.nrao.edu 




Fig. 1. — A composite of selected multi- wavelength observations 
of HLSW-Ol. Each image is 18" X 18". Clockwise from the up- 
per left; Keck Kp AO image; Spitzer/MIPS 24/xm; Spitzer/IKAC 
4.5/im; Subaru SuprimeCam i. In all panels the contours show the 
SMA 880/im observations. The numbering scheme used to identify 
sources in this letter is shown in the bottom left panel. The high- 
est resolution Herscfcei/SPIRE passband has a resolution of 18.6", 
larger than the entire field shown here. 



TABLE 1 

Photometry 



Wavelength [ftm] 


Flux Density 


Telescope/Detector 


0.48 (g) 


26.0 ±0.3 AtJy 


INT/WFC 


0.63 (r) 


46.5 ±0.6 /iJy 


INT/WFC 


0.76 (j) 


47.9 ± 0.3 /iJy 


Subaru / SuprimeCam 


2.2 [Kp) 


63.1 ± 2.0 AtJy 


KeckII/NIRC2 


4.5 


376 ± 6 /iJy 


Spitzer/IKAC 


5.8 


442 ± 11 ^tjy 


Spitzer/IKAC 


8.0 


558 ± 16 AtJy 


Spitzer/IKAC 


24 


5.5 ± 0.4 mJy 


Spitzer/MIPS 


72 


22.2 ± 3 mJy 


Spitzer/MIPS 


160 


310 ± 8 mJy 


Spitzer/MIPS 


250 


425 ± 10 mJy 


Herschel/SPIKE 


350 


340 ± 10 mJy 


Herschel/SPIKE 


510 


233 ± 11 mJy 


Herschel/SPIKE 


880 


52.8 ±0.5 mJy 


SMA 


1000-1100 


27.5 ±0.6 mJy 


CSO/Z-Spec 


1100-1200 


20.4 ±0.5 mJy 


CSO/Z-Spec 


1200-1300 


16.2 ±0.5 mJy 


CSO/Z-Spec 


1300-1400 


12.0 ±0.5 mJy 


CSO/Z-Spec 


1400-1500 


9.9 ± 0.6 mJy 


CSO/Z-Spec 


3400 


0.61 ±0.19 mJy 


CARMA 


214000 


1.8 ± 0.7 mJy 


VLA 



Note. — Combined flux densities for all the detected images 
as detailed in the text. For brevity, only the summed flux density 
from all images is provided. All values are calibrated relative to a 
= const SED. Note that calibration errors, which are dominant 
at most wavelengths and are strongly correlated between points, 
are not included, and neither are contamination/confusion errors. 
These values have not been corrected for magnification. 

not high enough resolution to explore this issue. There- 
fore, we assume that the near- and far-IR emission are 
co-located. Gil find evidence for a small shift (~ 0.4") 
between the PdBI CO emission and the optical emission, 
too small to significantly affect our results. 
The standard form for a modified blackbody SED 
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is a frequency-dependent optical depth factor: cx 
(1 — exp [— r (v)]) Bi, (T), where is the Planck func- 
tion. The optical depth is assumed to be a power 
law in frequency, r = {v/vq)^ following iDraind ()2006[ ). 
Ao = c/vq is the wavelength where the optical depth is 
unity. In the optically thin case, i/ <C I'd, this reduces to 
oc {T). The latter is often used in the literature 

to estimate temperatures, but here our data allow us to 
drop the assumption of optical thinness. In both cases 
we join the modified blackb ody to a s imple power law on 
the blue side of the peak (.Blain et al..,2003') . which only 
affects the 70/im observation. 

We fit both models from 70/^m to 1.5mm to de- 
rive the temperature and total IR luminosity, includ- 
ing the error on the magnification. The fit is shown 
in the left hand panels of Figure [21 and the parame- 
ters are given in Table [2] All of the parameters are 
well-constrained by our data. We find Aq ~ 200//m, 
a reasonable match to the theoretically expected value 
Ao — 100/xm (IDraind 120 06') and similar to that derived 
for Arp 220 (j Blain et a l. 2003). The emission is opti- 
cally thick bluer than observer-frame ^ 800/xm. A two- 
temperature model decreases the by < 0.002 for the 
optically thick model (since there is virtually no contri- 
bution from the second component for the best fit), and 
for the optically thin model by about 5, so the latter 
remains a very poor fit. 

The temperatures and (3 values for the two different 
models disagre e strongly; this is a lso the case for the fits 
to Arp 220 in iBlain et all ()2003[ ). The general model 
(i.e., 1 — exp [— r]) fits our data quite well, but the opti- 
cally thin model does not, with a reduced > 3. The 
derived temperature for the former is fairly high (~ 90 
K), suggesting a dust-enshrouded active galactic nucleus 
(AGN) contribution. This result is robust against remov- 
ing the data from any single instrument or any individual 
data point. Obtaining an acceptable fit for a more typ- 
ical dust temperature, such as 60K, requires increasing 
all the errors by a factor of > 2.2, including calibration 
errors. It is possible that the lensing is selectively magni- 
fying a warm component, and that this high temperature 
is not representative of the SMG as a whole. The poor 
quality of the optically thin fit is similarly robust, unless 
the SPIRE observations are removed, in which case it be- 
comes acceptable (x^ = l-l)- It is therefore possible that 
such warm SMGs have been missed in previous surveys 
that did not have observations near the peak of the SED, 
as it would then be difficult to distinguish between the 
optically-thin and thick cases. Since the general model 
fits the data much better, and makes fewer assumptions, 
henceforth we only discuss the results of this fit. 

The above findings are independent of the lensing 
model, unless the location of the emission (and hence the 
magnification) changes strongly from 250/xm to 1.5mm, 
which is unlikely. Turning to quantities which must 
be corrected for the lensing magnification, we find a 
far-IR luminosity of Lir = 1.43 x 10^'^ Lq, where Ljr 
is defined as the luminosity from 8 to 1000/im in the 
rest frame. This implies a star formation rate of ~ 
2500Moyr-i from the relation of lKennicutg (|1998[ ). ig- 
noring any AGN contribution. We also measured Ljr by 
spline-interpolating the observations, which gives a sim- 
ilar value (1.49 x 10^'^ L0). Assuming a mass-absorption 



TABLE 2 
Modified Blackbody Fits 



Model 


General: 1 — e 


Optically Thin: 


T 


88.0 ±2.9 K 


48.5 ±2.6 K 


P 


1.95 ± 0.14 


1.61 ±0.15 


Ao 


197 ± 19 ^ln^ 


NA 


LiR 


(1.43 ±0.09) X 10^^ Lq 


(1.13 ±0.09) X 10^^ Lq 


SFR 


2460 ± 16OM0yr-i 


1950 ± 160 M0yr-i 




6.77 for 6 dof 


22.1 for 7 dof 



Note. — Fit values for the two modified blackbody models, 
applied to the magnification corrected 70 — 1500/xm photometry. 
The second model assumes optically thin emission, and is only 
presented for comparison with literature values. The derived pa- 
rameters include the uncertainty in the magnification. For Lir, we 
assume h = 0.7, f2m = 0.27, and f2j\ = 0.73. 




1.0 1.5 2.0 2.5 3.0 1 10 100 1000 

p Wavelength [^m] 



Fig. 2.— SED fits to HLSW-01 after correction for flux mag- 
nification. The left-hand panels show the modified blackbody 
fits to the long-wavelength observations (top panel) and the con- 
straints on the temperature and /3 (the contours correspond to 
68/95/99% enclosed probability) for the general form of the model. 
Note that the data points are strongly correlated, and the opti- 
cally thick fit is an excellent match to the data. The right hand 
panel compares the optical to millimeter SED compared with the 
best fitting CIGALE model, the best combination of templates from 
Rowan- Robinson et al. (2010), and the IR SEDs of Arp220 and 
M82 (prototypical nearby IR-luminous galaxies) scaled to have the 
same maximum flux density. All wavelengths are observer frame. 

coefficient of = 2.64 m^kg^^ at 125/im following 
iDunne et al.l (|2003[ ). from the temperature and luminos- 
ity we infer a dust mass of Md ~ 1 x 10^ M©; k^, and 
hence M^, is uncertain by at least a factor of three. Fur- 
ther as suming a molecula. r gas to dust ratio of 60 for 
SMGs (jCoppin et al.ll200"8h . we estimate a gas depletion 
time of ~ 2.4 x 10^ yr, considerably shorter than the 
value of > 4 x lO'' yr for 'typical' SMGs derived by e.g., 
IDaddi et al.l (l20Tol) : see Rll for actual gas mass esti- 
mates. 

This source is weakly detected at 1.4 GHz in the FIRST 
survey, so we also compute qir, the logarithmic ratio of 
LiR and the rest-frame 1.4 GHz flux density. HLSW-01 
has an moderately low ratio, ^ir, = 1.5 ± 0.2, compared 
with the mean value and scatter of gjR = 2.4 ± 0.12 
for HerME S sources with firm radio cross-identifications 
(jivison et a l. 2010a), although not the lowest found. The 
high 1.4 GHz fiux density suggests some AGN contribu- 
tion. However, the 5" resolution of the FIRST survey is 
not good enough to rule out foreground contamination. 
While this is likely small, giR should be regarded as a 
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lower limit. 

In addition to the above simple models, we have also 
investigated a variety of template fits across all wave- 
lengths using several packages. A similar study was car- 
ried out for isolated , but somewhat int rinsically brighter, 
Herschel sources bv lBrisbin et al.l (j2010.' ). who found that 
model templates were generally a good match to the 
data. This is not the case here; available templates gen- 
erally do not fit both the near-IR and far-IR through 
millimeter-wave data simultaneously. The templates of 
[Rowan- Robinson et al. (2008, 2010) underestimate the 
fa r-IR flux by ^ 30%. This i s also the case for the models 
of iSiebenmorgen &: Kriigell (|2007D . In both cases this is 
because the SED of HLSW-01 pe aks blueward of the tem - 
plates. Similarly, the models of Ida Cunha et aTl (|2008[ ) 
are unable to reproduce the full SED (E. da Cunha 2010, 
private communication). The potential foreground con- 
tamination discussed in Section 2 is too small to explain 
these issues. 

Next we turn to the CIGALE package ()Noll et al.ll2009D . 
which combines optical/near- IR templates with a longer- 
wavelength dust model. Using a two-stellar population 
model, this fits the SED considerably better, but has 
some issues in the near-IR, overpredicting Kp (rest frame 
g) flux density by a factor of two and missing the slope 
of the IRAC observations. We are unable to explain this 
discrepancy. Ignoring these issues, CIGALE finds a total 
stellar mass of logj^p ^* — lO-^^o 3' ^-^d a star formation 
rate of log;^g SFR = S.Slgg, in M0 and M0 yr~^, respec- 
tively. 70 ± 30% of the stars are in a young, strongly 
extinguished stellar component {Ay = 6 ± 2), whose age 
is not well constrained. A sample of SED fits is shown 
in Figured] 



4. CONCLUSIONS 

The unusually large image separation of HLSW-Ol 
compared with most lensed sub-millimeter sources pro- 
vides an opportunity to study a sub-millimeter galaxy 
at a level of detail typically only possible for more lumi- 
nous sources. Due to confusion noise, this will generally 
not be feasible for un-lensed sources until the completion 
of ALMA. Detailed models of the gas and dust content 
based on CO emission are presented in Sll and Rll. 

Modified blackbody fits to the long-wavelength data 
{70fim to 1.5mm) imply a warm dust temperature of 
90 K, and a star formation rate of ^ 2500 M0 yr^^ . Com- 
pared with other SMGs, we find a short gas depletion 
timescale of only a few million years, assuming negligi- 
ble AGN contribution to Lir. This is one of the few 
SMGs that have been studied across such a wide wave- 
length range, due to the large image separation, so it is 
interesting that SED fits from the optical to the millime- 
ter are generally a fairly poor fit to the data, typically 
matching the short wavelength data well but underpre- 
dicting the far-IR peak. We obtain somewhat better re- 
sults with CIGALE, but it significantly overpredicts the 



2/im fiux and is not a great match to the IRAC obser- 
vations. It is unclear if this galaxy is simply unusual, or 
if the templates — which were largely derived from even 
brighter galaxies selected at longer wavelengths — are 
not a good representation of galaxies selected at wave- 
lengths near the peak of the far-infrared background. 
The latter would have significant implications for the in- 
ferred history of high-z star formation. Such warm sys- 
tems may have been missed by previous surveys lacking 
data near the peak of the SED because of the common 
assumption of optical thinness. Models suggest there 
should be a large number of strongly lensed SMGs in 
Hersc hel data, which is c onsistent with early observa- 
tions (jNegrello et al.|[20Tol) . so we expect to address this 
question soon. The variations in lensing magnification 
make this a promising tool to study SMGs across a wide 
range of intrinsic luminosities, although relatively few 
sources will have such a large image separation. Ad- 
ditional multi-wavelength observations, particularly at 
high resolution, would improve the SED constraints sig- 
nificantly. 

SPIRE has been developed by a consortium of insti- 
tutes led by Cardiff Univ. (UK) and including Univ. 
Lethbridge (Canada); NAOC (China); CEA, LAM 
(France); IFSI, Univ. Padua (Italy); lAC (Spain); Stock- 
holm Observatory (Sweden); Imperial College London, 
RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); Gal- 
tech, JPL, NHSC, Univ. Colorado (USA). This develop- 
ment has been supported by national funding agencies: 
CSA (Canada); NAOC (China); CEA, ONES, CNRS 
(France); ASI (Italy); MCINN (Spain); SNSB (Sweden); 
STFC (UK); and NASA (USA). The SubmiUimeter Ar- 
ray is a joint project between the Smithsonian Astro- 
physical Observatory and the Academia Sinica Institute 
of Astronomy and Astrophysics. The IRAM Plateau 
de Bure Interferometer is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). The Na- 
tional Radio Astronomy Observatory is a facility of 
the National Science Foundation operated by Associ- 
ated Universities, Inc. Support for CARMA construction 
was derived from the Gordon and Betty Moore Founda- 
tion, the Kenneth T. and Eileen L. Norris Foundation, 
the James S. McDonnell Foundation, the Associates of 
the California Institute of Technology, the University of 
Chicago, the states of California, Illinois, and Maryland, 
and the National Science Foundation. Ongoing CARMA 
development and operations are supported by by NSF 
grant ATL0838178 to CARMA, and by the CARMA 
partner universities. The authors would like to thank 
Elisabete da Cunha for running her models for us. The 
Herschel data presented in this letter will be released 
through the Herschel Database in Marseille, HeDalvF^. 

Facilities: Herschel (SPIRE), CSO (Z-Spec), Subaru 
(SuprimeCam), SMA, Hale (SWIFT), Keck:II (NIRC2), 
IRAM: Interferometer, ING: Newton (WFC), Spitzer 
(IRAC), Spitzer (MIPS), CARMA, VLA, GBT (Zpec- 
trometer) 
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